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ABSTRACT 


The  Synthetic  Bathymetric  Profiling  System  (SYNI3APS) 
consists  of  10  FORTRAN  IV  computer  programs,  a  random-access 
storage  device,  and  an  initial  bathymetric  data  base  of  over 
3  million  data  points.  SYNBAPS  is  designed  for  rapid 
generation  of  random  omnidirectional  bathymetric  profiles 
in  digital  form  along  great—circlfj  paths.  The  initial  data 
base  will  cover  most  of  the  Northern  Hemisphere  and  will  be 
extended  to  other  regions  as  suitable  bathyi?ietric  contour 
charts  become  available. 

Data  derived  from  the  bathymetric  contour  charts  are 
structured  into  a  griuded  data  surface  by  the  application  of 
a  cubic  spline  algorithm-  The  gridded  data  are  stored  on  a 
random-access  storage  device  by  5-degree-squ<ire  areas.  An 
accessing  program,  initiated  by  a  user's  request,  extracts  the 
5-degree-square  blocks  of  data  for  processing.  The  interpolation 
of  the  final  profile  is  accomplished  by  orienting  a  cubic  spline 
algorithm  along  a  great-circle  path  and  interpolating  the  depth 
values  from  the  5-degree  squares  falling  on  the  path.  A  status 
program  checks  the  content  and  condition  of  the  random-access 
storage  device. 

SYNBAPS  will  provide  bathymetric  profiles  at  about  one-fifth 
the  cost  and  one-hundredth  the  time  of  present  semiautomated 
methods. 
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FOREWORD 


This  report  describes  a  computer  system  and  programs  that 
will  establish  a  world-wide  bathymetric  data  bank  and  g=*nerato 
computer-drawn  bathymetric  profiles.  The  research  was  performed 
by  the  Naval  Oceanographic  Office  in  support  of  the  Office  of 
Naval  Research,  Long  Range  Acoustic  Propagation  Project,  wr. ich 
provided  funding.  It  is  part  of  a  major  bathymetric  char*  ir<'‘ 
project  covering  the  North  Atlantic  and  North  Pacific  Oceans. 
Bathymetric  data,  usually  in  the  form  of  profiles,  are  essntial 
elements  in  the  development  of  acoustic  propagation  models  '.nd 
predictions,  which  are  required  for  naval  planning,  systers 
development,  and  operations.  The  computerized  bathymetric 
profiling  system  and  specialized  data  bank  described  here  will 
generate  computer-drawn  bathymetric  profiles  at  a  small  fraction 
of  the  time  and  cost  of  manually  produced  profiles.  This 
specialized  data  bank  will  be  operational  when  approximately 
600  5-degree-square  areas  have  been  structured  on  a  random-accej.s 
storage  device.  Presently,  the  contour  data  required  for  the 
structuring  procedure  are  being  d; gitized  under  ONR-LRAPP 
contract  No.  N00014-72-C-0466. 


P.  V.  PURKRABEK 
Captain,  U.S.  Navy 
‘O'.runander 

U.S.  Naval  Oceanographic  Office 
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INTRODUCTION 


The  need  for  a  computerized  bathymetric  data  bank  and 
techniques  for  rapidly  manipulating  large  quantities  of  data 
became  evident  as  demand  upon  the  Naval  Oceanographic  Office 
for  bathym.etric  profiles  increased  and  became  more  urgent.  It 
became  increasingly  difficult  to  satisfy  these  demands  through 
manual  compilation  of  depth  soundings,  contouring,  and  profile 
constructions.  A  massive  recompilation  and  reanalysis  of 
bathymetric  data,  systematic  revision  of  all  bathyn>etric  charts 
in  the  North  Atlantic  and  North  Pacific  Oceans,  including 
extension  of  chart  coverage  to  the  equator,  was  underway.  At 
the  same  time  the  impracticality  of  using  the  existing  data  bank 
of  bathymetric  soundings  for  machine  generation  of  profiles 
became  apparent.  The  need  for  a  specialized  bathymetric  data 
bank  to  support  acoustic  -  oceanographic  modeling  gave  rise  to 
development  of  a  synthetic  bathymetric  profiling  project  using 
the  new  bathymetric  contour  charts  as  the  data  base.  The  project 
developed  procedures  for  digitizing  the  contour  charts,  and 
computer  programs  and  subroutines  for  data  storage  and  retrieval 
and  for  profile  generation.  Mr.  Thomas  M.  Davis,  Naval  Oceano¬ 
graphic  Office,  provided  special  assistance  in  developing 
programs  SPLINT  (SYNGRID) ,  BURNS  (SYNC0N2R) ,  BATHY  (subroutine 
BATHY)  and  DAV7HAT  (SYNCHEX)  and  contributed  to  the  basic  philosophy 
regarding  SYNBAPS.  Mr.  J.D.  Brown,  Naval  Oceanographic  Office, 
assisted  in  the  software  development  and  digitization  of  the  test 
data. 


Funds  for  this  project  were  provided  by  the  Office  of  Naval 
Research  through  the  Long  Range  Acoustic  Propagation  Project. 

One  of  the  basic  inputs  to  most  Navy  long-range,  acoustic 
propagation  models  are  bathymetric  profiles  in  digital  form.  These 
profiles  usually  are  plotted  along  a  great-circle  path  (glossary) 
as  a  function  of  range  versus  depth.  Two  methods  of  generating 
such  profiles  generally  have  been  employed.  In  the  first,  a  ship 
sails  a  predetermined  great-circle  path  collecting  continuous 
bathymetry  using  a  precision  depth  recorder  (PDR) .  if  the  course 
is  acvurately  adhered  to,  the  PDR  record  can  ne  merged  with  the 
navigational  record  to  obtain  the  bathymetric  profile.  If  the 
navigational  record  is  poor,  the  track  of  the  ship  will  have  to 
be  adjust' 1  and  normalized  to  obtain  a  satisfactory  bathymetric 
profile.  A  profile  thus  produced  is  accurate  and  retains  most 
of  the  high  frequency  information  but  is  costly  in  ship  time, 
hard  to  schedule,  and  usually  results  in  only  a  single  profile. 


A  second  means  of  obtaining  a  bathymetric  profile  is  to  plot 
a  great-circle  path  on  a  bathymetric  contour  chart,  or  series  of 
charts,  and  digitize  the  range  and  depth  at  the  intersection  of 
the  path  with  each  bathymetric  contour.  When  a  large  number  of 
great-circle  profiles,  each  several  thousand  miles  long,  involving 
dozens  of  bathymetric  charts,  are  constructed,  the  labor  costs  are 
considerable.  Profiles  produced  manually  from  charts  tend  to  be 
schematic,  blocky,  and  subject  to  human  error.  Most  importantly, 
both  of  these  methods  are  slow  and  cannot  be  achieved  in  real  time. 

Although  various  phases  of  both  methods  have  been  automated, 
within  the  Navy  and  elsewhere,  no  totally  satisfactory  solution 
has  been  achieved  to  the  present  time.  The  system  proposed  in 
this  report  is  one  approach  to  solving  the  above  problems. 

The  Synthetic  Bathymetric  Profiling  System  (SYNBAPS)  is  a 
combination  of  digital  computer  software  (programs)  and  a  random- 
access  storage  file  (presently  a  CDC  813  permanent  disk)  of  gridded 
bathymetric  data,  employed  to  generate  random,  great-circle, 
bathymetric  profiles  suitable  for  acoustic  propagation  modeling. 
SYNBAPS  is  completely  automatic,  requiring  only  the  input,  via  a 
control  card,  of  the  latitude  and  longitude  of  the  beginning  and 
end  points  to  extract  the  desired  profile.  The  profile  also  can 
be  generated  given  the  latittide  and  longitude  of  the  beginning 
point,  the  bearing,  and  the  maximum  range.  The  generated  profile 
is  available  in  two  forms.  The  first  is  a  computer-drawn  profile 
where  range  in  whole  nautical  miles  is  plotted  against  depth,  in 
either  meters  or  fathoms;  the  second  is  a  punched  card  deck  of 
the  same  data.  The  profile  outputs  in  card  image  are  available  on 
magnetic  tape  where  large  quantities  of  data  are  involved. 

A  bathymetric  profile  along  a  great-circle  path  of  about 
8,000  nm  can  be  generated  in  approximately  3  minutes  of  computer 
time  on  a  second  generation  computer  and  can  be  plotted  in  about 
3  minutes  on  an  incremental  plotter.  A  cost  comparison  shows  that, 
by  present  semiautomatic  methods,  a  set  of  19  short  profiles 
totaling  9,000  nm  required  144  man-hours  at  a  cost  of  $900.  The 
same  profiles  could  be  produced  by  SYNBAPS  in  1.4  man-hours  at  a 
cost  of  $50,  for  a  savings  of  18;1  in  dollars  and  100:1  in  time. 

OUTLINE  OF  SYSTEM  OPERATION 

The  SYNBAPS  software  can  be  broken  down  into  three  distinct 
program  functions  associated  v/ith  structuring,  accessing,  and 
status  (fig.  1) . 

The  structuring  programs  create  a  gridded  bathymetric  data 
base  and  structure  it  on  a  random-access  device  in  precise 
form.  The  smallest  cell  of  the  data  base  is  a  5-minute-square 
grid  where  the  north-south  side  is  in  meridional  minutes  or  parts 
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FIGURE  1.  SYNTHETIC  BATHYMETRIC  PROFILING  SYSTEM  DIAGRAM 
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and  the  east-west  side  is  in  longitudinal  minutes.  On  a 
Mercator  projection  contour  chart  this  is  a  5-ininute  rectangular 
grid.  The  bathym»jtric  data  are  logically  formatted  to  place 
depth  values  at  the  intersection  of  each  5-minute  grid  crossing 
as  shown  in  figure  2. 

The  next  level  of  structuring  is  to  index  the  5-minute  cells 
into  5-degree  squares  called  Marsden  Square  Locator  numbers  (MSQLOC) 
using  the  Marsden  square  system  which  divides  the  earth  surface 
into  5-degree  squares  (fig.  3) .  Further  subdivision  of  the 
Marsden  square  by  quadrants  is  shown  in  figure  4.  The  MSQLOC  is 
the  quadrant  number  followed  by  the  Marsden  square  number  as 
follows: 


Marsden  square  number+quadrant  =  MSQLOC 
Example:  036+2  =  0362 

The  MSQLOC  is  a  unique  v/orldwide  reference  to  each  5-degree  square 
of  gridded  bathymetric  data.  The  MSQLOC  area  includes  a  5-minute 
overlap  of  all  sides  as  shown  in  figure  5  for  MSQLOC  0362. 

The  gridded  bathymetric  data  base  is  created  following  the 
procedure  used  by  Davis  and  Kontis  (1970) .  However,  accurate 
synthetic  data  derived  from  large  and  medium-scale  bathiTTietric 
charts  are  used  instead  of  original  survey  data.  The  synthetic 
track  data  are  derived  from  charts  by  superimposing  parallel  track 
lines,  5  minutes  apart,  over  the  MSQLOC  area.  Extraction  of  the 
data  usually  starts  from  the  lower  left  corner.  The  orientation 
of  the  track  lines  can  be  any  direction  from  v/est-east  (90°  bearing) 
to  nearly  south-north  (1°  bearing) ,  but  not  true  north,  which 
necessitates  changing  several  statements  in  the  gridding  program. 
The  only  other  restriction  is  that  the  first  track  be  a  west-east 
track  across  the  MSQLOC  area.  The  remaining  tracks  may  be  of  any 
orientation  and  in  any  order. 

The  data  are  extracted  from  the  chart  by  digitising  the 
intersections  of  the  synthetic  track  v;itn  the  contours  sequentially 
along  the  track.  Interpolated  points  must  be  extracted  for  the 
beginning  and  end  of  each  full  track.  These  tracks  must  extend 
5  minutes  beyond  the  MSQLOC  area  on  all  sides  as  shown  in  figure  6. 
Short  tracks  may  be  added  to  emphasize  certain  topographic 
characteristics  such  as  spot  elevations.  These  can  be  extracted 
at  any  orientation  except  true  north-south  as  shown  in  figure  6B. 

Each  digitized  track  is  assigned  a  sequence  number,  but  the 
physical  order  of  the  tracks  in  the  card  deck  is  arbitrary  after 
the  first  track.  These  digitized  tracks  are  inputs  to  the  gridding 
program.  The  output  from  that  program  is  a  punched  deck  of  gridded 
bathymetric  data  with  the  point  or  origin  in  the  lower  left  corner. 
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FIGURE  3.  MARSDEN  SQUARE  CHART 
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FIGURE  5.  EXAMPLE  OF  MSQLOC  AREA 
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EXAMPLE  OF  SYNTHETIC  TRACK  ORIENTATION 


These  data  are  physically  unformatted.  A  number  of  error 
checks  are  made  before  and  after  the  gridded  bathymetric  data 
are  created.  The  gridded  data  then  are  placed  on  a  random-access 
storage  device  using  a  predetermined  "look-up"  table  (list  of 
acronyms) .  At  this  point  the  data  are  ready  to  be  accessed. 

At  present,  a  bathymetric  profile  can  be  generated  up  to 
8,000  nm  long  and  crossing  30  MSQLOC  areas.  This  limitation  can 
be  increased  if  necessary.  The  accessing  is  initiated  by 
supplying  the  latitude  and  longitude  of  a  beginning  and  end 
point  or  the  latitude  and  longitude  of  a  beginning  point  with  the 
bearing  and  maximum  range.  Combinations  of  these  accessing  schemes 
can  also  be  used. 

The  first  step  in  retrieving  a  profile  from  the  data  bank 
is  to  generate  its  great-circle  path.  At  the  same  time  each 
MSQLOC  area  that  the  path  crosses  is  identified  and  a  search 
table  of  MSQLOC  areas  is  created.  For  each  MSQLOC  area  the 
search  table  contains  the  latitude  and  longitude  of  the  first  and 
last  point  in  that  MSQLOC  area,  the  forward-looking  bearing  at 
both  points,  the  accumulated  range  from  zero  for  both  points, 
and  the  MSQLOC  area  number.  In  turn,  each  MSQLOC  area  is  called 
from  the  random-access  storage  device  via  the  "look-up"  table 
and  the  profile  for  that  block  of  data  is  generated.  This  partial 
profile  is  then  placed  on  a  temporary  magnetic  tape.  The  next 
MSQLOC  area  is  called  from  the  random-access  storage  device  and 
the  cycle  is  repeated.  At  the  end  of  profile  generation  the 
temporary  magnetic  tape  is  rewound.  The  plotting  program  is 
then  called  and  the  partial  profiles  are  linked,  punched  on  cards, 
and  plotted  and/or  written  on  magnetic  tape. 

The  accessing  program  is  structured  so  that  long  profiles 
generally  are  processed  faster  than  numerous  short  profiles  that 
total  the  same  mileage.  The  great-circle  path  generation  requires 
about  20  seconds  for  an  8,000-nm  profile  plus  about  5  seconds 
for  each  full  MSQLOC  area  crossed  for  the  interpolation. 

The  only  maintenance  to  be  performed  to  the  system  is  the 
eventual  updating  of  the  gridded  bathymetric  data  based  on  the 
random-access  storage  device.  This  is  easily  accomplished  by 
recycling  through  the  structuring  phase  of  the  system  any  MSQLOC 
area  that  requires  updating  and  then  replacing  that  block  on  the 
random-access  storage  device. 

A  status  report  can  be  generated  to  check  any  or  all  MSQLOC 
areas.  This  report  includes  the  random-access  device's  compatible 
data  block  size,  the  actual  column  and  row  sizes,  the  date  the 
data  block  was  added  to  the  random-access  device,  the  MSQLOC 
area  number,  the  relative  address,  and  the  actual  data,  if  required. 


10 


SOURCE  MATERIAL 


Bathymetric  contour  charts  instead  of  recorded  water  depths, 
are  the  source  for  the  SYNBAPS  data  base.  No  computer  algorithm 
(glossary)  that  can  successfully  handle  all  qualities  of  bathy- 
metric-track-line  data,  resolve  all  navigational  errors,  and  can 
apply  a  contouring  philosophy  to  such  data  has  been  developed. 
These  functions  require  the  subjective  judgement,  based  on 
knowledge  of  geologic  processes,  of  the  bathymetrist  whose  final 
product  is  the  bathymetric  contour  chart.  The  bathymetrist ' s 
very  subjectivity  creates  the  data  continuity  which  is  a  requisite 
element  of  SYNBAPS.  A  long  profile  requires  on  omnidirectional, 
continuous  data  base,  something  that  is  seldom  achieved  with 
either  survey  or  random  ship  track  line  data  alone.  Using  areas 
having  high  quality  and  dense  data  coverage  as  a  framework,  the 
bathymetrist  extends,  interpolates,  and  extrapolates  regional 
trends  into  areas  of  lesser  data  to  build  a  continuous  picture 
of  the  submarine  topography. 

Although  SYNBAPS  is  designed  for  worldwide  application, 
initially  a  data  base  will  be  created  only  for  the  Northern 
Hemisphere,  and  possibly  the  Indian  Ocean.  Other  regions  will 
be  added  to  the  data  base  when  sufficient  continuous  data  become 
available.  The  charts  used  for  the  North  Pacific  Ocean  will  be 
large  to  medium  scale  (1:1,000,000  or  larger)  versions  of  the 
U.S.  Naval  Oceanographic  Office  H.O.  Pubs.  1301,  1302,  and  1303 
(U.S.  Naval  Oceanographic  Office  1969,  1971A  and  B) .  Recent 
unpublished  large-to-medium  scale  charts  compiled  by  the  U.S. 

Naval  Oceanographic  Office  will  be  utilized  for  the  North  Atlantic 
and  Mediterranean  Sea.  Where  applicable,  classified  data  can  be 
incorporated  in  the  data  base  without  compromising  security.  The 
gridded  data  point  from  a  classified  chart,  which  was  contoured 
from  classified  data  or  from  a  mixture  of  classified  and  unclassi¬ 
fied  data,  V7ill  be  indistinguishable  from  a  data  point  from  an 
unclassit’ed  chart.  Only  the  originator  will  know  which  depth 
values  were  created  from  classified  data  and  that  they  may  be 
more  accurate  than  other  points.  The  originator  will  keep  a 
separate  noncomputerized  file,  indexed  by  MSQLOC  areas,  showing 
the  source  of  the  contours,  their  evaluation,  classification, 
and  other  pertinent  information.  There  will  be  no  reference  to 
original  track  spacing,  area  limits,  navigation,  sounding  device, 
or  platform  within  the  data  base.  Preparation  of  the  charts  for 
digitization  is  discussed  in  more  detail  in  appendix  A. 

SYSTEM  DESCRIPTION 

A.  Structuring  Programs 
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The  relationship  between  structuring  programs  is  given  in 
a  flow  diagram  in  figure  7.  The  main  processing  programs  are 
SYNTRACK,  SYNCARD,  SYNCHEX,  SY^1GRID,  SYNC0N2R,  and  SYNBLOCK  (list 
of  acronyms) .  One  additional  program  that  is  unique  to  this 
pa.ticular  system  is  the  digitizer  scaling  program  (CALMA  485) 
which  scales  on  a  Mercator  chart  the  latitude,  longitude,  and 
depth  for  each  contour  intersection  along  the  track.  The  output 
from  this  program  is  a  binary  magnetic  tape  of  scaled  values. 

Any  digitizer  and/or  digitizer  processor  program  can  be  used  as 
long  as  it  generates  the  same  program  elements  regardless  of 
output  mode. 

The  MSQLOC  area  to  be  digitized  is  mounted  and  leveled  on 
the  digitizer  table  (fig.  7)  .  Starting  in  the  lovirer  left  corner 
each  track  is  scanned  for  data  points  from  left  to  right  and 
from  bottom  to  top.  The  tracks  are  scanned  an  additional  5 
minutes  on  each  end  to  permit  interpolation  rather  than  extra¬ 
polation  on  end  points  in  the  gridding  program.  The  MSQLOC 
identification  and  operator  name  are  entered  as  a  header  informa¬ 
tion  group  before  the  u?ta  scanning  is  begun.  The  binary  coded 
decimal  (BCD)  magnetic  i-ape  generated  by  the  digitizer  is 
processed  by  the  CALMA  4f-5  processor  program  to  produce  a  binary 
magnetic  tape  of  scaled  latitude,  longitude,  and  depth  data.  The 
binary  tape  is  processed  by  SYNTRACK  which: 

•  breaks  up  the  data  string  into  tracks, 

•  checks  for  missing  data  points, 

•  checks  for  operator  errors, 

•  reformats  the  data  to  card  image,  and 

•  punches  out  a  header  card,  track  card,  data  cards 
(one  point  per  card) ,  and  a  blank  card. 

An  illustration  of  this  deck  structure  is  given  in  figure  8.  After 
errors  have  been  corrected,  the  card  deck  generated  by  SYNTRACK  is 
run  through  the  SYNCARD  program.  This  program  checks  to  insure 
that  the  longitudes  of  contour  intersections  are  not  repeated,  but 
either  increase  or  decrease  depending  upon  quadrant.  In  addition, 
this  program  tests  the  depth  value  to  determine  if  it  is  within 
about  plus  or  minus  two  times  the  contoiir  interval.  In  regions 
of  rapid  depth  change  contours  may  be  skipped  if  they  are  evenly 
spaced.  All  errors  are  flagged  for  correction. 

After  all  corrections  have  been  made,  the  card  deck  is  run 
through  the  track  plotting  program  (SYNCHEX) .  This  program  plots 
the  tracks  as  they  were  digitized  and  annotates  each  contour 
intersection  on  the  synthetic  track  line  with  cross  ticks.  This  plot 
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FIGURE  7.  SYNBAPS  STRUCTURING  PROGRAMS  FLOW  DIAGRAM 


ETC. 


FIGURE  8.  OUTPUT  DECK  STRUCTURE  FROM  SYNTRACK 
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insures  that  the  proper  and  sufficient  number  of  points  have 
been  extracted  from  the  MSQLOC  area.  Additional  tracks  of  data 
can  be  created  at  this  time,  if  required  by  the  complexity  of 
the  submarine  topography.  SYNCHEX  requires  a  control  card  that 
is  in  reality  the  first  data  card.  The  format  for  this  card  is 
given  in  figure  9.  If  no  further  corrections  or  additions  are 
to  be  made  to  the  data  deck,  the  MSQLOC  area  is  ready  for 
conversion  to  gridded  bathymetry. 

The  SYNGRID  program  is  fundamental  to  the  structuring  phase 
of  SYNBAPS.  SYNGRID  transforms  the  synthetic  track  line  data 
into  gridded  bathymetric  data.  The  mathematical  foundation  and 
philosophy  behind  the  one-dimensional  cubic  spline  used  to 
structure  the  gridded  data  is  fully  explained  by' Davis  and  Kontis 
(1970) .  SYNGRID  is  a  modification  by  Davis  of  his  original 
program  (SPLINT)  to  handle  bathymetric  data  instead  of  gravity 
data.  SYNGRID  is  very  flexible  as  it  grids  track-line-point 
data  on  either  a  Mercator  projection  or  a  Cartesian  coordinate 
system  and  can  compute  mean  data  for  various  size  cells  on  either 
system.  Summarizing  Davis  and  Kontis  (1970; ,  the  value  of  this 
method  lies  not  only  in  its  ability  to  fit  the  observed  data 
values  but  to  retain  the  continuity  of  the  first  and  second 
derivatives.  This  method  might  be  considered  the  mathematical 
analog  of  the  draftsman's  plastic  spline. 

Because  the  cubic  spline  is  a  function  of  only  one  independent 
variable,  the  data  obtained  along  a  synthetic  track  line  must  be 
adjusted  to  lie  on  a  straight  line.  Under  most  conditions  this 
creates  no  problem  as  the  data  are  digitized  along  straight  lines. 
The  interpolation  formula  used  by  Davis  fits  each  data  exactly, 
has  continuous  first  and  second  derivations,  and  is  a  simple  cubic 
polynomial  in  x  within  the  interval  between  each  pair  of  data 
points.  The  distance  along  the  track  line  then  may  be  interpreted 
as  the  independent  variable.  Therefore,  taking  the  data  from  one 
track  at  a  time,  the  position  of  the  data  points  are  converted 
into  X,  y  coordinates  and  a  least  squares  straight  line  is  fitted 
to  these  locations.  Because  no  statistical  significance  is 
attached  to  this  operation,  either  x  or  y  may  be  considered  the 
independent  variable.  The  computer  program  listed  in  appendix  B 
considers  x  the  equivalent  longitude  as  the  independent  variable. 

If  the  survey  tracks  happen  to  run  exactly  north-south,  the 
program  should  be  modified  to  consider  y  as  the  independent 
variable. 

The  perpendicular  distance  be  ’een  the  least  squares  straight 
line  and  each  data  point  is  determined  and  used  to  project  the 
points  orthogonally  onto  the  line  with  an  adjusted  data  value 
(based  on  an  estimate  of  the  local  gradient)  assumed  to  be  a 
function  of  distance  only.  If  the  perpendicular  distance  between 
this  point  and  the  least  squares  line  is  less  than  predetermined 
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NLINE  XMIN  XMAX  YMIN  YMAX  XiN  YIN 


I  II  21  31  41  51  61  71 

1  FORMAT  (I  10,  6  F  10.0) 

R=  RIGHT  JUSTIFIED  •  ^FLOATING  POINT  REQUIRED 


NLINE  =  TOTAL  NUMBER  OF  TRACKS 

XMIN  =  minimum  number  of  minutes  from  prime  meridian  -  east  or  west 

XMAX  =  maximum  number  of  minutes  from  prime  meridian  -  east  or  west 

YMIN  =  minimum  number  of  meridional  ports  from  equator  -  north  or  south 

YMAX  =  maximum  number  of  meridional  parts  from  equator  -  north  or  south 

XIN  =  east*wesr  dimension  of  plot  in  inches 
YIN  =  north-south  dimension  of  plot  in  inches 

NOTE:  1 .  North  and  east  are  positive,  south  and  west  are  negative. 

2.  MSQLOC  areas  require  5°  minutes  of  overlap  on  all  sides. 

3.  Meridional  ports  are  found  in  reference; 

Naval  Oceanographic  Office,  1962 
H.O.  Pub.  9 -Table  5 


FIGURE  9.  SYNCHEX  CONTROL  CARD  FOR  TRACK  PLOTTING  OF  MSQLOC 
AREA 


pivot  distance  (usually  set  at  0.2  of  a  meridional  part),  the 
value  associated  with  the  data  point  is  unchanged.  If  this 
distance  is  greater  than  the  pivot  distance,  then  the  adjusted 
value  associated  v;ith  the  mapped  coordinates  is  computed. 

In  the  computer  program  for  the  cubic  spline  algorithm 
(SPLINE)  contained  in  appendix  B,  the  pivot  distance  is  selectable 
via  a  control  card.  This  pivot  distance  is  usually  equal  to 
the  maximum  distance  which  one  could  move  a  data  point  v;ithout 
significantly  changing  its  value.  In  order  to  minimize  the  error 
associated  with  the  assumption  that  the  gradient  correction  is 
independent  of  direction,  continuous  synthetic  survey  tracks 
which  deviate  appreciably  from  a  stright  line  should  be  broken 
up  into  smaller  segments  v?ith  each  segment  treated  as  a  separate 
track.  The  mapped  coordinates  and  adjusted  data  values  may 
be  considered  as  irregularly  spaced  digital  samples  from  a 
function  whose  independent  variable  is  distance  along  the  track 
from  some  arbitrary  starting  point,  and  whose  dependent  variable 
is  the  adjusted  data  values. 

Utilizing  the  mapped  data,  the  cubic  spline  is  determined 
for  each  track.  The  cubic  spline  may  then  be  used  to  interpolate 
data  values  at  the  intersections  of  the  straight  least  square 
track  lines  and  a  set  of  parallel  lines  whose  spacing  is  equal 
to  the  desired  final  grid  spacing  (5  minutes) .  If  the  direction 
of  the  survey  tracks  is  predominantely  east-west  then  the  direction 
of  the  set  of  parallel  lines  is  north-south.  Similarly,  for 
north-south  tracks,  the  lines  are  run  east-west. 

The  computer  program  (app.  B) ,  is  designed  to  operate  on 
tracks  in  any  direction,  except  exactly  north-south.  The 
direction  of  the  set  of  parallel  grid  lines  is  controlled  by 
the  direction  of  track  line  number  one.  Since  the  track  number 
designation  is  arbitrary,  this  feature  allows  the  user  to 
determine  the  desired  orientation  (W-S  or  E-K)  of  the  parallel 
grid  lines  in  order  to  obtain  as  many  intersections  as  possible. 

The  interpolated  data  values  generated  as  outlined  in  the 
preceding  paragraph  may  be  regarded  as  unequally  spaced  digital 
samples  from  a  function  whose  independent  variable  is  distance 
along  each  of  the  parallel  lines.  Application  of  the  spline 
procedure  in  this  cross  track  direction  produces  the  final 
interpolated  values  at  the  desired  grid  points.  If  mean  anomalies 
are  desired,  grid  points  are  generated  at  one-half  the  final 
grid  spacing  and  the  resulting  nine  points  are  averaged  to 
produce  the  mean  value  for  each  grid  cell. 

The  control  card  formats  for  SYNGRID  are  given  in  figure  10. 
The  output  from  SYNGRID  is  a  new  punched  card  deck  of  gridded 
bathymetry  with  seven  points  per  card.  The  printout  from  SYNGRID 
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ISET 


ISETS  =  number  of  MSOLOC  areas  fo  be  processed  during  a  computer  run 


MEAN 


ALAT  =  latitude  of  MSQLOC  for  lower  lefthand  comer  in  degrees 

ALONG  =  longitude  of  MSQLOC  for  lower  lefthand  comer  in  degrees 

PLAT  =  latitude  of  MSQLOC  for  upper  righthand  comer  in  degrees 

PLONG  =  longitude  of  MSQLOC  foi  upper  righlhand  comer  in  degrees 
GRID  =  grid  spacing  for  output  dota  in  minutes 

MEAN  =  blank,  no  mean  computed;  =1,  mean  computed 

ITOT  =  total  number  of  trocks  of  input  data 

IT YPE  =  1 ,  grid  is  in  Mercator  prelection;  =  "I ,  grid  is  in  X  ond  Y  units 

PIVOT  =  maximum  distance  from  track  for  pivot  test 

MSQLOC  =  Marsden  Square  Locator  orea  number 


FIGURE  10.  SYNGRID  CONTROL  CARDS  FOR  GRIDDING  TRACK  DATA 
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V’i.ll  indicate  if  the  MSQLOC  area  has  been  structured  correctly. 

An  even  more  efficier.c  method  of  chejking  is  to  pass  the  gridded 
bathymetric  data  through  the  syNC0N2R  program. 

The  SyNC0N2R  program  (fig.  7)  plots  contours  of  the  gridded 
data  on  a  Mercator  projection  at  the  same  scale  as  the  source 
manuscript.  The  source  manuscript  can  be  overlain  by  the 
gridded-data  contour  plot,  for  a  comparison  of  content  and 
form.  This  plotting  check  requires  a  29-inch  drum  plotter  or 
equivalent,  while  the  SYNC0N2R  program  itself  requires  a  control 
card  (fig.  11) .  In  addition,  the  DATA  statement  variable  (CL) 
requires  a  specification  of  the  contour  levels  that  will  be 
plotted  (see  app.  B) .  An  optional  DATA  statement  variable 
(LABELS)  can  be  used  if  labels  are  desired  (see  app.  B) .  If  the 
SYNC0N2R  plot  is  satisfactory,  the  gridded  bathymetry  is  loaded 
on  the  random-access  storage  device  via  the  loading  program 
(SYNBLOCK),  fig.  7). 

Before  a  block  of  gridded  data  can  be  loaded  on  the  random- 
access  storage  device,  the  device  must  be  primed  with  a  traffic 
director  program  (SYNTABLE,  fig.  7) .  SYNTABLE  is  a  predetermined 
"look-up"  +-able,  which  gives  SYNBLOCK  basic  information  that  is 
needed  to  place  a  block  of  gridded  data  in  its  proper  address  on 
the  device.  Using  the  MSQLOC  area  number  as  the  key,  the  table 
supplies  the  relative  address,  the  actual  block  size  to  be 
transmitted,  and  a  file  key  or  name.  The  file  key  indicates  by 
name  in  which  file  in  the  storage  device  a  particular  block  of 
data  is  to  be  placed.  ?.n  example  of  the  "look-up"  table  printout 
is  given  in  table  1.  In  the  DATA  statement  N  is  equal  to  the 
number  of  MSQLOC  areas  now  on  the  "look-up"  table.  The  relative 
address  is  the  physical  location  from  the  beginning  of  the  file 
of  the  first  v/ord  of  the  data  block.  The  actual  block  size  is  the 
quantity  of  storage  required  to  contain  the  data  plus  the 
identification  groups  and  is  an  even  multiple  of  32  (Aiken,  et  al. 
1970).  The  storage  requirement  for  the  actual  block  size  is 
predetermined  and  is  listed  in  table  2  by  hemisphere  latitude 
bands,  which  include  the  overlap. 

Using  the  "look-up"  table  from  SYNTABLE  on  the  random-access 
storage  device,  a  block  of  gridded  bathymetry  can  nov/  be  loaded 
by  SYNBLOCK.  The  punched  deck  of  gridded  data  is  preceded  by 
two  header  cards.  The  first  card  contains  the  number  of  sets  to 
be  loaded  and  the  second  card,  one  for  each  set,  specifies  the 
MSQLOC  area  number  and  the  column  and  row  information  obtained 
from  table  2  (see  app.  B  for  exact  card  formats) .  The  DATA 
statement  N  is  equal  to  the  number  of  MSQLOC 's  presently  on  the 
"look-up"  table.  SYNBLOCK  then  looks  up  the  file,  the  relative 
address,  and  the  block-size  information  from  the  preloaded  table 
for  each  MSQLOC  area  and  places  the  data  in  its  proper  location. 

An  identification  group  containing  the  following  is  placed  at  the 
end  of  the  data  block: 
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I  5  9  13  17  21 

1  FORMAT  (5  I  4,  4  F  10.0) 
— ^  R  =  RIGHT  JUSTIFIED 


• 

• 

• 

31 

41 

51 

•  =  FLOATING  POINT  REQUIRED 


N 

= 

number  of  columns  of  Input  doto 

M 

s= 

number  of  rows  of  Input  doto 

NCL 

= 

number  of  contour  levels 

MM 

s 

row’s  maximum  array  size 

NN 

= 

column's  maximum  array  size 

XA 

= 

1.0  =  minimum  number  of  rows 

YA 

= 

1.0  =  minimum  number  of  columns 

XG 

s 

X  axis  width  of  plot  In  Inches 

YG 

s 

y  axis  height  of  plot  In  Inches 

FIGURE  11 .  SYNCON2R  CONTROL  CARD  FOR  CONTOUR  PLOTTING 
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SYNBAP5  DISK  FILE  LOCATOR  TABLE 


MSQLOC 

RELATIVE 

ADDRESS 

SIZE  OF 
BLOCK 

211 

0 

3936 

212 

3936 

3936 

213 

7672 

4000 

214 

11872 

4000 

571 

15872 

4064 

572 

19936 

4064 

573 

24000 

4128 

574 

28128 

4128 

931 

32768 

4256 

932 

37024 

4256 

933 

41280 

4448 

934 

45728 

4448 

1291 

50176 

4704 

1292 

54880 

4704 

1293 

59584 

4928 

1294 

65536 

4928 

1651 

70464 

4928 

1652 

75776 

5312 

1653 

81088 

5824 

1654 

86912 

5824 

2011 

98304 

6464 

2012 

104768 

6464 

2013 

112096 

7328 

2014 

119424 

7328 

TABLE  I.  EXAMPLE  OF  "LOOK  UP"  TABLE  ROM 


FILE 

KEY 


E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 

E08C 


SYNTABLE 
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TABLE  II.  SYNBAPSWORD  STORAGE  REQUIREMENTS 


NUM  =  actual  size  of  storage  block 
ICOL  =  number  of  columns  of  array 

IROW  =  number  of  rows  or  array 

MSQLOC  =  Marsden  Square  Locator  area  number 
IDAY  =  day  that  data  ware  placed  in  storage 
MONTH  =  month  that  data  were  placed  in  storage 
lYEAR  =  year  that  data  were  placed  in  storage 
LOCATE  =  relative  address 

This  completes  the  structuring  phase  of  SYNEAPS.  The 
punched  cards  of  gridded  bathymetric  data  are  loaded  on  magnetic 
tape  with  one  MSQLOC  area  per  file  using  a  UTILITY  program 
(Rozanski,  et  al.  1968) .  This  magnetic  tape  is  saved  for  backup 
to  the  random-access  storage  device. 

B.  Accessing  Programs 

The  relationship  between  accessing  programs  is  given  in  a  flow 
diagram  in  figure  12.  The  tvra  accessing  programs  are  SYNBAPSl  and 
SYNPLOT  (app.  C)  .  The  rtBquest,  in  the  form  of  control  cards,  is 
submitted  to  the  SYNBAPSl  program  (fig.  13) .  The  formats  for  this 
request  may  be  either  all  "BEARINGS"  or  all  "POINTS"  or  can  be 
a  mixture  of  both,  as  long  as  the  number  of  beams  is  correctly 
indicated  for  each  set  (the  variable  NOOFBM) . 

With  the  exception  of  SYNGRID,  only  a  brief  explanation  of 
the  program's  operation  was  given  in  the  structuring  phase 
discussion.  Because  SYNBAPSl  and  SYNPLOT  may  be  used  by  others, 
they  will  be  described  in  more  detail. 

Figure  14  contains  a  more  detailed  program  flov/  diagram  of 
SYNBAPSl.  VJhen  a  request  is  submitted  to  SYNBAPSl  the  first 
operation  is  to  call  in  the  SEAARCH  subroutine  to  generate  the 
great-circle  path  to  be  followed  by  the  profile.  SE/vARCH  uses 
both  the  direction  solution  of  the  great  circle,  subroutine 
GCDIST,  and  the  indirect  solution  GCPATH  (Chang,  1969A  and  B)  to 
create  a  latitude,  longitude,  forward  bearing,  and  range  for 
each  nautical-mile  point  from  the  beginning  to  the  end  of  a 
profile.  In  addition,  subroutine  MSQFQ  is  used  to  calculate  the 
MSQLOC  area  for  each  of  the  points.  SEAARCH  then  creates  a  range 
search  table  of  only  those  points  that  start  a  profile,  enter 
or  exit  a  MSQLOC  area,  or  terminate  a  profile.  This  table  is 
printed  out  and  also  placed  in  COMMON. 
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PRINT  OUT 
OF  RESULTS 
OF  PROFILING 


FIGURE  12.  SYNBAPS  ACCESSING  PROGRAMS  FLOW  DIAGRAM 


24 


■f  / 

f;, 


1  FORMAT  (  •  5,  A8) 

NOOFBM  =  number  of  profiles  of  NCARD  type  to  be  processed 
NCARD  =  "BEARINGS"  OR  "POINTS" 

L^-=LEFT  JUSTIFIED  -►R=RIGHT  JUSTIFIED 


IDBM  ALAfl'^AN^fiLMINY  BS 


I  9  12  16  18  21  2527  Z7  47 

1  FORMAT  (A6,  2X,  2((F3.0,  F3.0),  IX,  A1,  IX),  2  F  10.0) 

FOR  NCARD  =  "BEARING" 

L<t-=LEFT  JUSTIFIED  -►Rs  RIGHT  JUSTIFIED  •  sFLDATING  POINT  REQUIRED 


AMIN  ALONG  AE  BMIN  BLONG  BE 
ID  ALAT  1  AN  UlMINIBLATI  BN  IbLMIN 


Sisi  I  jsss  i  ISB!  I  SIS!  I 


9  12  16  18  21  2527  30  3436  59  4345 

2  FORMAT  (A6,  2X,  4(F3.0,  F3.0,  IX,  Al,  IX)) 


FOR  NCARD  =  "POINTS" 


IDBM  =  unique  profile  i.D.  (alphanumeric) 

ALAT  =  degree  of  lotitude  -  start  point 

AMIN  =  minute  of  latitude  *  stort  point 

AN,  BN  -■  hemisphere  indicator,  N  or  S 

ALONG  =  degree  of  longitude  -  start  point 

ALMIN  =  minute  of  longitude  -  stort  point 

AE,  BE  =  hemisphere  indicotor,  E  or  W 

BS  =  beoring  from  start  point  for  "BEARING"  cord  only 

DD  =  maximum  ronge  from  stort  point  for  "BEARING"  cord  only 

BLAT  =  degree  of  latitude  -  end  point 

BMIN  =  minute  of  latitude  -  end  point 

BLONG  =  degree  of  longitude  -  errd  point 

BLMIN  =  minute  of  longitude  “  end  point 


FIGURE  13.  SYNBAPS1  PROFILE  REQUEST  CONTROL  CARD 
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FIGURE  14.  ACCESSING  PROGRAMS  DETAIL  FLOW  DIAGRAM 
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Subroutine  MINCON  is  called  in  to  calculate  the  starting 
point  for  the  profile  within  MSQLOC  in  minutes  from  the  lower 
left  corner.  MIMCON  uses  the  function  AMP  to  calculate  the 
meridional  parts  for  the  latitude  component.  The  mathematical 
foundation  for  AMP  is  given  in  Thomas  (1964)  and  in  U.S.  Naval 
Oceanographic  Office  (1962) . 

Subroutine  RHUMB  is  called  in  to  calculate,  using  AMP,  the 
rhumb  line  bearing  through  the  MSQLOC  area.  A  rhumb  line  is 
used  here  because  the  subroutine  BATHY  can  only  interpolate 
along  a  straight  line.  The  rhumb  line  approximates  a  chord  of 
the  great-circle  path  cn  a  Mercator  chart  with  the  maximum  deviation 
from  the  great  circle  at  the  approximate  midpoint  of  that  chord 
in  the  MSQLOC  area.  This  deviation  varies  from  zero  to  a 
maximum  of  about  two  nautical  miles  depending  upon  the  great- 
circle  path  orientation.  Maximum  deviations  occur  in  east-west 
paths  in  high  latitudes,  but  are  considered  a  necessary  trade¬ 
off  for  the  system's  overall  speed  of  operation. 

The  random-access  storage  device  is  queried  by  the  subroutine 
LOOKUP,  v;hich  passes  through  the  SYNTABLE  to  find  the  file  key 
and  the  relative  address  of  the  MSQLOC  area,  then  extracts  the 
actual  block  size  and  the  column  and  row  information.  These 
parameters  are  used  by  the  subroutine  BATHY  to  extract  gridded 
bathymetric  data  for  the  MSQLOC  area. 

From  subroutine  BATHY  the  subroutine  GRIDBLK  calls  in  the 
gridded  data.  Subroutine  BATHY  determines  which  quadrant  the 
rhumb  line  will  pass  through  so  as  to  maximize  the  number  of 
intersections  for  interpolation.  This  quadrant  will  determine 
whether  or  not  the  columns  or  the  rows  will  be  the  independent 
variable  for  tr ''  cubic  spline.  The  quadrant  arrangement  is 
shown  in  figure  15. 

If  the  rhumb  line  falls  in  quadrants  2  or  4 ,  the  direction 
of  the  first  interpolation  is  along  a  column  and  the  independent 
variable  is  the  distance  from  the  origin  along  the  column  to  the 
intersection  of  the  rhumb  line.  If  the  rhumb  line  falls  in 
quadrants  1  or  3,  the  interpolation  will  be  along  a  rov/  and  the 
independent  variable  then  is  the  distance  from  the  origin  along 
the  row  to  the  intersection  with  the  rhumb  line.  At  the  inter¬ 
section  a  value  is  interpolated  by  the  cubic  spline  using  the 
gridded  data  values  along  that  column  (or  rov;)  as  the  dependent 
variable. 

When  all  the  values  have  been  interpolated  at  each  inter¬ 
section,  the  values  now  become  the  dependent  variable  v/hile  the 
distance  along  the  rhumb  line  from  the  start  of  the  profile 
becomes  the  independent  variable.  The  cubic  spline  is  used  once 
more  to  interpolate  the  final  profile  values  at  distances 
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MSQLOC  ORIGIN 


FIGURE  15.  QUADRANTS  FOR  SUBROUTINE  BATHY 
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corresponding  to  every  meridional  part  along  the  rhumb  line  to 
the  end  of  the  MSQLOC  area.  An  example  of  this  rotation  is 
given  in  figure  16.  When  a  profile  for  a  MSQLOC  has  been 
generated,  BATHY  calls  the  PUNOUT  subroutine  to  put  the  MSQLOC 
profile  data  on  a  temporary  magnetic  tape.  MERPIX  and  AMP  are 
used  by  PUNOUT  to  calculate  the  rhumb  line  distance  in  meridional 
parts  and  set  up  a  scaling  factor.  The  parameters  are  used  by 
PUNOUT  to  adjust  the  profile  generated  by  BATHY,  which  is  in 
meridional  parts  versus  depth,  to  a  profile  which  shows  nautical 
miles  versus  depth  by  linear  interpolation.  Only  when  these 
operations  are  complete  is  the  MSQLOC  profile  data  v/ritten  on 
the  temporary  magnetic  tape  and  the  next  MSQLOC  area  or  the  next 
profile  processed. 

Each  segment  of  a  profile  represents  a  single  MSQLOC  area. 

When  the  individual  segments  are  written  on  the  temporary  magnetic 
tape  the  depth  is  in  the  same  units  as  in  the  gridded  data  base 
and  the  range  is  in  'nautical  miles  from  starting  point  within 
the  MSQLOC  area,  which  in  each  case  is  zero.  At  the  end  of  the 
SYNBAPSl  program  the  temporary  magnetic  tape  is  rewound.  The 
program  SYNPLOT  then  reads  this  tape  either  on  the  same  or  a 
subsequent  run.  As  each  MSQLOC  profile  segment  is  read  into 
SYNPLOT  it  is  linked  in  sequence  to  the  other  MSQLOC  areas  to 
produce  a  great-circle  profile.  If  geometric  conversion  to  other 
depth  units  is  required,  it  is  performed  at  this  point. 

When  the  great-circle  profile  is  complete,  it  is  punched 
out  on  cards  and  the  profile  is  plotted.  This  process  is 
repeated  for  as  many  profiles  as  desired.  Although  the  format 
for  the  punched  profile  cards  is  fi.xed  at  eight  depth-versus- 
range  points  per  card,  the  profile-plotting  format  is  very 
flexible.  This  flexibility  is  attained  through  a  control  card 
for  SYNPLOT,  the  format  for  which  is  given  in  figure  17.  Generally, 
whenever  SYNBAPSl  cannot  find  a  MSQLOC  block  of  gridded  data  on 
the  random-access  storage  device  or  the  plotting  dimensions  are 
not  set  for  minimal  limits  (fig.  17) ,  the  processing  will  halt 
at  that  point  and  skip  to  the  next  profile,  allowing  the  job  run 
to  continue  while  an  error  message  is  printed  out. 

The  profiles  generated  by  SYNBAPS  are  intended  as  input  to 
long-range,  acoustic  propagation  models.  Although  not  necessarily 
accurate  to  geophysical  or  geodetic  standards,  the  sythetic 
profiles  are  interpolated  to  the  accuracy  required  by  the  models. 

A  depth  value  is  interpolated  at  each  nautical -mile  point  from 
the  starting  point  to  the  terminus  of  the  profile  along  a 
great-circle  path.  Latitude  and  longitude  values  are  rounded 
to  the  nearest  minute,  and  the  range  is  rounded  to  the  nearest 
nautical  mile. 
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FIGURE  16.  PROFILE  EXTRACTION  ROM  GRIDDED  DATA  BASE 


I  UNIT 


YLTH  CONVERT 


e 

• 

L4- 

□ 

• 

• 

1 

II 

21 

28 

31 

41 

51 

100  FORMAT  (  2  F  10.0,  A  7,  3  X,  2  F  10.0) 

L4-  =  LEFT  JUSTIFIED  •  =FL0ATIN6  POINT  REQUIRED 


R 

D 

lUNITS 

YLTH 


CONVERT 


I 


X  axis  scaling  factor,  nautical  miles  per  inch 
y  axis  scaling  factor,  meters  or  fathoms  per  inch 
lobel  for  y  axis  (x  axis  is  always  in  nautical  miles) 
the  total  height  of  the  y  axis  plot  that  will  be  displayed,  the 
maximum  is  10  inches.  This  usually  set  as  o  multiple  of  D. 
Example:  when  plotting  at  500  fathoms/ inch,  to ‘be  able  to 
display  a  profile  that  goes  down  to  a  depth  of  4500  fathoms, 

YLTH  would  equal  9  inches.  If  not  correctly  set  or  if  plot 
exceeds  13  feet  on  the  x  axis,  that  profile  Is  omitted  from 
plotting.  However,  the  cards  are  still  punched, 
the  data  base  is  uncorrected  for  speed  of  sound  in  sea  water. 

For  fathoms  the  assumed  standard  is  800  fathoms  per  second,  for 
meters  it  is  1500  meters  per  second.  To  convert  from  fathoms 
to  meters  CONVERT  =  1.8750,  from  meters  to  fothoms  CONVERT 
0.533“”"3.  If  no  conversion  needed  CONVERT  =  blank  or  0.0. 


FIGURE  17.  SYNPLOT  CONTROL  CARD 
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The  great-circle  subroutines  are  based  upon  a  sphere 
21,600  nm  in  diameter  and  can  have  a  maximum  error  of  20  nm 
over  a  distance  of  1  hemisphere  (about  11,000  nm) .  This 
amounts  to  an  error  of  about  2nm/l,000nm  of  range.  For 
profiles  of  1,000  nm  or  less  this  error  is  insignificant  in 
propagation  model  applications,  but  it  could  be  important  at 
very  long  ranges.  The  magnitude  of  this  error  depends  upon 
the  difference  in  shape  between  the  sphere  and  the  oblate 
spheroid  and  on  the  method  of  path  generation.  Greater  accuracy 
can  be  obtained  by  using  a  geodesic  v/here  the  error  is  1  m  in 
latitude,  longitude,  and  range  and  0.035  sec.  in  bearing  within 
3  hemisphere  (Thomas,  1965  and  1970) . 

Within  each  MSQLOC  area  there  is  a  difference  between  the 
path  follov^ed  by  the  great  circle  and  the  actual  path  along  which 
the  depths  values  are  interpolated  (fig.  18)  .  Because  SYNBAPSl 
requires  a  straight  line  along  which  to  interpolate  depth  values, 
a  rhumb  line  between  the  first  position  entering  a  5-degree  square 
and  the  last  position  before  leaving  the  square  is  used  instead 
of  the  curved  great-circle  path.  For  all  great  circles  that 
follow  a  meridian  or  the  equator  this  difference  is  zero.  For 
all  other  directions,  the  maximum  difference  is  located  at  the 
approximate  mid -point  along  a  rhumb  line  within  5-degree  square. 
Under  the  most  unfavorable  condition  of  high  latitude  and  an 
east-west  orientation,  this  difference  rarely  exceeds  2  nm. 

Preliminary  estimates  of  the  accuracy  of  the  interpolated 
depth  values  in  the  profile  plane  are  +  15  fm.  This  assumes  that 
there  are  no  positional  errors  in  the  great-circle  path  in  the 
horizontal  plane.  A  completed  data  bank,  including  regions  of 
smooth  to  rough  topography,  will  be  needed  before  full  error 
analysis  can  be  undertaken. 

C.  Status  Program 


Program  SYNSTAT  queries  the  random-access  storage  device 
through  the  SYdTABLE  for  a  listing  of  the  identification  group 
from  each  .MSQLOC  gridded  daca  block.  This  listing  includes  the 
file  key  as  in  the  following  example; 

ACTUAL  DATE  ADDED  TO 


MSQLOC 

FILE 

KEY 

RELATIVE 

ADDRESS 

BLOCK 

SIZE 

NO.  OF 
COLUMNS 

NO.  OF 
ROWS 

RANDOM-ACCESS 

DEVICE 

1. 

1291 

E08C 

50176 

4704 

63 

74 

18  April  1972 

2. 

1292 

E08C 

54880 

4704 

63 

74 

19  April  1972 

All  MSQLOC  gridded  data  blocks  or  selected  ones  can  be  listed. 
They  are  selectable  through  SYNSTAT  control  cards  as  shown  in 
figure  19. 
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LATITUDE 


LONGITUDE 


FIGURE  18.  DlFF£REh:CE  BETWEEN  RHUMB  LINE  AND  GREAT 
CIRCLE  PATH  WITHIN  A  FiVE-DEGREE  SQUARE 
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ITYPE  NUM 


mm 

a 

- ^ 

1 

8 

15 

10  FORMAT  (A7,  I  7) 

1  L<-  =  LEFT  JUSTIFIED  — ►R  =  RI6HT  JUSTIFIED  ] 

iTYPE  ~  "ALL",  the  contents  cf  the  complete  random  occess  storage  device 
will  be  listed. 

"PARTIAL",  only  those  MSQLOC's  listed  on  the  following  control  cords 
will  be  listed  out. 

NUM  =  if  blank,  all  the  MSQLOC’s  listed;  if  present  only  that  number  of  MSQLOC's 
on  the  following  control  cords  will  be  listed. 


IA(1) 


1A(I6) 


A 


H 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

-♦R 

6 

It 

16 

21 

26 

31 

36 

41 

46 

Oi 

56 

61 

66 

71 

76 

I 

20  FORMAT  (16  I  5) 
ALL  RIGHT  JUSTIFIED 


R  =  RIGHT  JUSTIFIED 


lA  =  array  of  MSQLOC  numbers. 


FIGURE  19.  SYNSTAT  CONTROL  CAROS 


D.  CDC  3800  System  Subroutines  and  Functions 


The  subroutines  used  to  open  the  file,  position,  read,  and 
v;rite  on  the  CDC  813  permanent  disk  are  on-line  COMPASS  language 
routines  provided  by  the  Naval  Research  Laboratory,  Research 
Computation  Center  Staff  (Aiken,  et  al.,  1970).  These  subroutines 
are  DKOPEN,  DKLOCATE,  DKREAD,  and  DKWRITE.  The  subroutine  DATA 
is  an  off-line  COMPASS  language  routine  that  retrieves  the  integer 
day,  month,  and  year  from  the  computer's  internal  clock  (Houston, 
1969) .  The  function  TIMELEFT  is  an  on-line  COMPASS  language 
routine  that  retrieves  time  marks  from  the  computer's  internal 
clock.  It  is  used  to  time  various  phases  of  the  structuring  and 
accessing  programs  operation  (Shannon,  1968) . 

The  "jn-line  plotting  subroutines  PLOTS,  PLOT,  LINE,  SYMBOL, 
and  AXIS  are  FORTRAN  language  routines.  With  the  possible 
exception  of  LINE  and  AXIS  these  routines  are  part  of  the  standard 
Calcomp  plotter  package  (Gossett,  et  al.,  pending). 

Most  o"  the  previously  mentioned  subroutines  and  functions 
are  unique  to  the  NRL  CDC  3800  computer  system.  However,  these 
routines  have  counterparts  on  any  large  computer  system,  and. 
their  replacement  should  pose  little  or  no  problem. 

PROFILE  OUTPUT 

Two  adjoining  MSQLOC  areas,  1291  and  1292,  in  the  western 
North  Pacific  Ocean  were  selected  to  test  the  computer  program 
and  vjere  digitized,  structured,  and  placed  on  the  random-access 
storage  device.  The  location  of  five  test  profiles  along  rhumb 
lines,  subsequently  shown  in  figures  21,  22,  and  23,  are  indexed 
in  figure  20.  The  contour  chart  used  as  an  index  chart  shows 
only  part  of  the  contour  data  that  will  input  to  the  data  base; 
therefore,  the  test  profiles  show  a  slight  difference  in  detail. 
Figure  2i,  a  profile  through  both  MSQLOC  areas,  shows  that  the 
link  point  between  two  data  blocks  is  undetectable.  This  530-nm 
profile  was  generated  in  7  seconds. 

In  figure  22,  composed  of  three  profiles  A,  B,  and  C,  a 
dashed  line  is  superimposed  on  each  profile.  The  dashed  lines 
are  profiles  hand  drawn  by  a  bathymetrist ,  and  the  solid  lines 
are  the  computer  profiles.  All  the  profiles  used  the  same  data 
base.  Although  the  general  shapes  for  both  types  of  profiles  are 
the  same,  the  cubic  spline  profiles  show  details  betv;een  the 
contour  levels  that  would  otherv;ise  be  lost  if  XK>t  captured  by  the 
surface  of  gridded  bathymetric  data.  This  is  especially  true  in 
the  more  steeply  sloping  areas  because  the  cubic  spline  considers 
data  adjoining  the  profile  path.  The  three  profiles  in  figure  22 
show  the  system's  ability  to  start  a  profile  inside  a  MSQLOC  area. 
Figures  22A  and  22C  show  profiles  that  terminate  in  gently  sloping 
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FIGURE  20.  INDEX  OF  SAMPLE  PROFILES 


irj 


areas,  and  figure  22B  shows  a  profile  terminating  on  the  upslope 
side  of  a  seamount  in  the  next  area  to  the  north.  Figures  22B 
and  22C  shov;  that  the  cubic  spline  can  follow  both  convex  and 
concave  submarine  topography  equally  well. 

Figure  23  shows  two  mirror-image  prof iles,  which  illustrate 
the  profile  repeatability  along  the  same  path  in  either  direction. 
Profile  A  was  run  from  west  to  east,  then  profile  B  was  run  from 
east  to  west,  both  along  the  same  path. 

FURTHER  MODIFICATIONS,  ADDITIONS  AND  OTHER  APPLICATIONS 

The  first  modification  to  SYNBAPS  v;ill  replace  the  great- 
circle  subroutines,  GCPATH  and  GCDIST,  in  the  accessing  phase 
with  geodesic  subroutines,  GEODIST  and  GEOPATH.  The  argument 
list  for  the  new  routines  will  be  the  same  as  for  the  great 
circle  routines.  The  second  modification  v;ill  replace  the  contour 
checking  program,  SYNC0N2R,  in  the  structuring  phase  with  a 
smoother  contour  plotting  program.  A  third  modification  will 
attempt  to  increase  the  overall  efficiency  (speed  of  operation) 
by  simplifying  the  programs.  One  example  is  to  use  buffering 
statements  when  writing  and  reading  on  the  temporary  magnetic 
tape. 


An  additional  program  to  operate  on  the  SYNBAPS  output  will 
be  an  updated  automatic  depth  correction  routine  based  on 
Matthews'  sound  velocity  correction  tables.  This  will  permit 
the  use  of  depth  values  either  corrected  or  uncorrected  for 
speed  of  sound  in  sea  water. 

An  additional  version  of  SYNBAPSl,  the  accessing  program, 
called  SYNBAPS2  is  be’ng  considered.  This  program  will  generate 
eight  radial  profiles  simultaneously  from  one  point  to  the  edge 
of  a  MSQLOC  area  or  an  irregular  chart  area.  This  output  could 
be  useful  for  profile  evaluation  of  site  locations  where  greater 
detail  is  required.  In  addition,  SYNBAPSl  can  be  merged  with  the 
NODC  Ocean  Station  Data  file  to  produce  a  composite  plot  of  the 
bottom  profile  and  selected  sound  velocity  profiles  along  a 
great-circle  path.  Extending  this  concept  one  more  step  will 
produce  profile  plots  of  various  acoustic  environmental  parameters, 
such  as  depth  to  the  axis  or  bottom  of  the  deep  sound  channel,  by 
marrying  SYNBAPS  to  an  appropriate  oceanographic  data  file  or 
files.  The  number  of  possibl.e  combinations  of  oceanographic  data 
with  the  depth  data  using  SYNBAPS  is  almost  infinite. 

A  system  similar  to  SYNBAPS,  but  using  land  topography, 
could  be  applied  in  radar  terrain  studies  and  weather  pattern 
models  requiring  elevation  data. 
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SUMMARY  AND  CONCLUSIONS 

The  SYNBAPS  data  base  was  designed  to  meet  the  specific 
and  immediate  need  for  bathymetric  profiles  for  acoustic 
modeling.  However/  properly  used/  it  offers  many  applications 
beyond  its  preliminary  designs. 

Often  in  naval  planning  as  well  as  in  naval  operations, 
speed  is  as  important  as  accuracy  when  information  is  needed. 
SYNBAPS  is  not  ideally  suited  to  hydrographic  charting  because 
some  high-frequency  information  is  lost,  but  it  provides  very 
rapid  responses.  SYNBAPS  has  these  additional  features: 

•  Only  data  points  are  stored  in  the  data  bank, 

•  The  locations  of  data  points  are  logically  structured 
on  a  Mercator  projection  by  5-minute  intersections, 

•  Random  access  to  the  data  is  by  large  blocks  (5- 
degree  square) , 

•  The  data  bank  is  updated  by  replacing  blocks  of  data, 

•  The  size  of  the  data  bank  is  fixed  once  it  has  been 
created  for  any  ocean  area, 

•  Classified  survey  data,  in  chart  form,  can  be  incor¬ 
porated  in  the  data  base  with  no  compromise  of 
security, 

•  Highly  compacted  forms  of  the  accessing  program  and 
the  data  bank  can  be  used  on  shipboard. 
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GLOSSARY  OF  SELECTED  TERMS 


Accuracy 

Address 

Algorithm 

Alphanumeric 

Argument  list 

Band  (Latitudinal 

Bathymetric 
Bathymetric  chart 


The  degree  of  freedom  from  error,  that 
is,  the  degree  of  conformity  to  truth 
or  to  a  rule.  Accuracy  is  contrasted 
with  precision,  e.g. ,  four-place 
numerals  are  less  precise  than 
six-place  numerals;  nevertheless  a 
properly  computed  four-place  numeral 
might  be  more  accurate  than  an 
improperly  computed  six-place  numeral. 

(1)  An  identification,  as  represented 
by  a  name,  label,  or  number,  for  a 
register,  location  in  storagv. ,  or  any 
other  data  source  or  destination  such 
as  the  location  of  a  station  in  a 
communication  network.  (2)  Loosely, 
any  part  of  an  instruction  that 
specifies  the  location  of  an  operand 
for  the  instruction. 

A  finite  set  of  rules  that  gives  a 
sequence  of  operations  for  solving  a 
specific  type  of  problem.  It  should 
have  the  following  features, 

(1)  Finiteness,  (2)  Definiteness, 

(3)  Input,  (4)  Output,  and 
(5)  Effectiveness. 

Pertaining  to  a  character  set  that 
contains  both  leti.ors  and  numerals, 
and  usually  other  characters. 
Synonymous  with  Alphameric. 

List  of  the  formal  parameters  of  a 
subprogram  used  as  an  explicit 
transfer  of  information  to  or  from  a 
subprogram. 

land)  Any  latitudinal  strip,  designated  by 
accepted  units  of  linear  or  angular 
measurement,  v;hich  circumscribes 
the  earth. 

Relating  to  the  measurement  of 
ocean  depths. 

A  topographic  map  of  the  floor  of 
the  ocean. 
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Bathymetry 


The  science  of  determining  and 
interpreting  ocean  depths  and 
topography . 


Bearing 


Binary 


Binary  Coded  Decimal 
(BCD) 


Block 


Block  diagram 


1.  (general)  The  horizontal  angle  at 
a  given  point  measured  clockwise 
from  a  specific  reference  datum  to  a 
second  point.  Also  called  bearing 
angle. 

2.  (navigational)  The  horizontal 
direction  of  one  terrestrial  point 
from  another,  expressed  as  the 
angular  distance  from  a  reference 
direction.  It  is  usually  measured 
from  000"  at  the  reference  direction 
clockwise  through  360°.  The  terms, 
bearing  and  azimuch  are  sometimes  used 
interchangeably,  but  in  navigation 
the  former  customarily  applies  to 
terrestrial  objects  an„  the  latter  to 
the  direction  of  a  point  on  the 
celestial  sphere  from  a  point  on  the 
earth 

(1)  Pertaining  to  a  characteristic  or 
property  involving  a  selection, 
choice,  or  condition  in  which  there 
are  two  possibilities.  (2)  per¬ 
taining  to  the  numeration  system  with 
a  radix  of  two. 

Pertaining  to  a  decimal  notation  in 
which  the  individual  decimal  digits 
are  each  represented  by  a  group  of 
binary  digits,  e.g.,  in  the  8-4-2-1 
binary  coded  decimal  notation,  the 
number  23  is  represented  as  0010  0011, 
whereas  in  binary  notation,  23  is 
represented  as  10111. 

A  set  of  things,  such  as  words, 
characters,  or  digits,  handled  as 
a  unit. 

A  diagram  of  a  system,  instrument, 
computer,  or  program  in  which  selected 
portions  are  represented  by  annotated 
boxes  and  interconnecting  lines. 
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Cartesian  coordinates 


COMMON 


COMPASS 

DAT/. 

Data 

Deck 

Dependent 

Digitize 


Values  representing  the  location  of  a 
point  in  a  plane  in  relation  to  two 
intersecting  straight  lines,  called 
axes.  The  point  is  located  by 
measuring  its  distance  from  each  axis 
along  a  parallel  to  the  other  axis. 

If  the  axes  are  perpendicular  to  each 
other,  the  coordinates  are  rectangular; 
if  not  perpendicular,  they  are  oblique 
coordinates.  This  system  is  extended 
to  represent  the  location  of  points  in 
three-  dimensional  space  by  referencing 
to  three  mutually  perpendicular 
coordinate  axes  which  intersect  at  a 
common  point  of  origin. 

Is  a  specification  statement,  used 
during  compilation  rather  than 
execution  as  a  convenient  method  for 
passing  values  between  main  program 
and  subprograms  without  mentioning 
them  as  arguments. 

Control  Data  Corporation  assembly 
language  for  CDC  3000-  and  6000- 
series  computers . 

Is  a  specification  statement,  used 
during  compilation  rather  than 
execution  as  a  convenient  method  for 
entering  data  value  into  referenced 
storage  areas. 

Any  representations  such  as  characters 
or  analog  quantities  to  v/hich  meaning 
might  be  assigned. 

A  collection  of  punched  cards. 

variable  A  fixed  variable  given  as  a  function 

of  another  variable,  i.e.,  if  y  is 
given  as  a  function  of  x,  then,  y  is 
the  dei>endent  variable. 

(1)  The  conversion  of  graphical  analog 
information  or  characters  into  digital 
form,  usually  for  the  purpose  of  rapid 
manipulation  or  storage  by  a  digital 
computer  (2)  to  express  data  in  a 
digital  form. 
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Field  In  a  record,  a  specified  area  used 

for  a  particular  category  of  data, 
e.g.,  a  group  of  card  columns  used  to 
represent  a  wage  rate  or  a  set  of  bit 
locations  in  a  computer  word  used  to 
express  the  address  of  the  operand. 

File  A  collection  of  related  records 

treated  as  a  unit.  Thus  in  inventory 
control,  one  line  of  an  invoice  forms 
an  item,  a  complete  invoice  forms  a 
record,  and  the  complete  set  of  such 
records  forms  a  file. 

Fixed  point  Pertaining  to  a  numeration  system  in 

which  the  position  of  the  point  is 
fixed  with  respect  to  one  end  cf  the 
numerals,  according  to  some  convention. 

Floating  point  Pertaining  to  a  numeration  system  in 

which  the  position  of  the  point  does 
not  remain  fixed  with  respect  to  one 
end  of  the  numerals. 

Flowchart  A  graphical  representation  for  the 

definition,  analysis,  or  solution  of  a 
probJem,  in  which  symbols  are  used  to 
represent  operations,  data,  flow  and 
equipment . 

Geodesic  A  line  of  shortest  distance  between 

any  two  points  on  any  mathematically 
defined  surface.  A  geodesic  line  is  a 
line  of  double  curvature,  and  usually 
lies  between  the  two  normal  section 
lines  which  the  two  points  determine. 
If  the  two  terminal  points  are  in 
nearly  the  same  latitude,  the  geodesic 
line  may  cross  one  of  the  normal 
section  lines.  It  should  be  noted 
that,  except  along  the  equator  and 
along  the  meridians,  the  geodesic  line 
is  not  a  plane  curve  and  cannot  be 
sighted  over  directly.  However,  for 
conventional  triangulation  the  lengths 
and  directions  of  geodesic  lines 
differ  inappreciably  from  corre¬ 
sponding  pairs  of  normal  section 
lines.  Also  called  geodesic  line; 
geodetic  line. 
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Great  circle 


A  circle  on  the  surface  of  the  earth, 
the  plane  of  which  passes  through  the 
center  of  the  earth. 

Header  card  The  first  card  or  cards  of  a  deck  of 

punched  cards  containing  identifi¬ 
cation  of  fixed  information  about  the 
punched  cards  of  variable  data  that 
follow. 

Independent  variable  A  variable  whose  assigned  value (s)  are 

arbitrary  when  defined  as  a  function 
of  another  variable,  i.e.  ,  if  y  is 
given  as  a  function  of  x,  then,  x  is 
the  independent  variable. 

Input  (1)  The  data  to  be  processed.  (2)  The 

stage  or  sequence  of  states  occurring 
on  a  specified  input  channel.  (3)  The 
device  or  collective  set  of  devices 
used  for  bringing  data  into  another 
device.  (4)  A  channel  for  impressing 
a  state  on  a  device  or  logic  element. 

(5)  The  process  of  transferring  data 
from  an  external  storage  to  an 
internal  storage. 

Interpolation  To  determine  intermediate  values 

between  given  fixed  values.  As 
applied  to  logical  contouring  to 
interpolate  is  to  ratio  vertical 
distances  between  given  spot  ."'levations. 

Lock-up  table  An  index  file  or  array (s)  which  is 

usually  used  to  access  a  main  record 
file.  It  contains  the  identifier  (or 
file  key)  and  the  storage  address  in 
sequential  or  non- sequential  order. 

It  may  also  contain  critical  infor¬ 
mation. 

MARSDEN  chart  A  system  introduced  by  Marsden  early 

in  the  nineteenth  century  for  showing 
the  distribution  of  meteorological 
data  on  a  chart;  especially  over  the 
oceans.  A  Mercator  map  proiection  is 
used;  the  world  between  90 '’I'  and  80 ®S 
being  divided  into  Marsden  "squares" 
each  of  10®  latitude  by  10®  longitude. 
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MARSDEN  chart  (Con.) 


Mercator  projection 


Merge 


Meridional  part 


Offline 


Online 


These  squares  are  systematically 
numbered  to  indicate  position.  Each 
square  may  be  divided  into  quarter 
squares,  or  into  100  1°  subsquares 
numbered  from  00  to  99  to  give  the 
position  to  the  nearest  degree. 

A  conformal  map  projection  of  the 
cylindrical  type.  The  equator  is 
represented  by  a  straight  line  true  to 
scale;  the  geographic  meridians  are 
represented  by  parallel  straight 
lines  perpendicular  to  the  line 
representing  the  equator;  they  are 
spaced  according  to  their  distance 
apart  at  the  equator.  The  geographic 
parallels  are  represented  by  a  second 
system  of  straight  lines  perpendicular 
to  the  family  of  lines  representing 
the  meridians  and  therefore  parallel 
with  the  equator.  Conf ormability  is 
achieved  by  mathematical  analysis, 
the  spacing  of  the  parallels  being 
increased  with  increasing  distance 
from  the  equator  to  confo.rm  with  the 
expanding  scale  along  the  parallels 
resulting  from  the  meridians  being 
represented  by  parallel  lines.  Also 
called  equatorial  cylindrical 
orthomorphic  map  projection. 

To  combine  two  or  more  sets  of  items 
into  one,  usually  in  a  specified 
sequence . 

The  length  of  the  arc  of  a  meridian 
between  the  equator  and  a  given 
parallel  on  a  Mercator  chart,  ex¬ 
pressed  in  units  of  one  minute  of 
longitude  at  the  equator. 

Pertaining  to  equipment  of  devices  not 
under  direct  control  of  the  central 
processing  unit. 

Pertaining  to  equipment  or  devices 
under  direct  control  of  the  central 
processing  unit. 
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Kandom  access 


The  dacreo  or  oi3cri~iriat.ior: 
v/hich  a  quanoioy  is  soared,  e.g.,  a 
ohree-digio  nun.erai  discriminates 
among  i,000  possibilities. 

A  veroical  secoion  of  one  surface  of 
the  ground,  or  rf  underlying  strata, 
or  both,  along  any  fixed  line. 

The  smallest  field  (group)  of  unique 
contiguous  characoers  or  digits. 

(1)  A  card  punched  v/ith  a  pattern  of 
holes  to  represent  data.  (2)  A  card 
as  in  (1)  before  being  punched. 

A  quantity  v/hose  successive 
integral  pov/ers  are  the  implicit 
irultipliers  of  the  sequence  of  digits 
that  represent  a  number.  For  example 
if  the  radix  is  5,  then  143.2  means 
1  times  5  to  the  second  power,  plus 
4  times  5  to  the  first  power,  plus  3 
times  5  to  the  zero  power,  plus  2 
times  5  to  the  minus  one  power. 

(1)  Pertaining  to  the  process  of 
obtaining  data  from,  or  placing  data 
into,  storage  where  the  time  required 
for  such  access  is  independent  of  the 
location  of  the  data  most  recently 
obtained  or  placed  in  storage.  (2) 
Pertaining  to  a  storage  device  in 
which  the  access  time  is  effectively 
independent  of  the  location  of  the 
data. 


Real  Time 


Relative  address 


Rhumb  line 


Round  off 


Routine 


Selection  overlay 


Storage 


(1)  Pertaining  to  the  actual  time 
during  which  a  physical  process 
transpires.  (2)  Pertaining  to  the 
performance  of  a  computation  during 
the  actual  time  that  the  related 
physical  process  transpires  in  order 
that  results  of  the  computation  can  be 
used  in  guiding  the  physical  process. 

Identifies  a  word  in  a  subroutine  or 
array  with  respect  to  its  position. 
Relative  addresses  are  translated 
into  absolute  addresses  by  the 
addition  of  some  specific  reference 
address,  usually  that  at  which  the 
first  word  of  the  routine  or  array 
is  stored. 

A  line  of  the  surface  of  the  earth 
making  the  same  angle  with  all 
meridians;  a  loxodrome  or  loxodromic 
curve  spiraling  toward  the  poles  in  a 
constant  true  direction.  Parallels 
and  meridians,  which  also  maintain 
constant  true  directions,  may  be 
considered  special  cases  of  the  rhumb 
line.  A  rhumb  line  is  a  straight  line 
on  a  Mercator  projection.  Also  called 
equiangular  spiral;  loxodrome,  loxo- 
droiaic  curve;  Mercator  track. 

To  delete  the  least  significant  digit 
or  digits  of  a  numeral  and  to  adjust 
the  part  retained  in  accordance  with 
some  rule. 

A  set  of  instructions  arranged  in 
proper  sequence  to  cause  a  computer  to 
perform  a  desired  task. 

A  tracing  of  selected  map  source 
detail  compiled  on  transparent 
material;  usually  described  by  the 
name  of  the  features  or  details 
depicted,  such  as  contour  overlay, 
vegetation  overlay.  Also  called  lift; 
pull  up;  trace. 

(1)  Pertaining  to  a  device  into  which 
data  can  be  entered,  in  which  it  can 
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W  V 


be  held,  and  from  which  it  can  be 
retrieved  at  a  later  time.  (2) 
Loosely,  any  device  that  can  store 
data.  (3)  Synonymous  v;ith  Memory. 

Produced  artifically;  devised, 
arranged,  or  fabricated  for  special 
situations  to  imitate  or  replace 
usual  realities. 
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LIST  OF  ACRONYMS  USED  IN  COMPUTER  PROGRAMS 


AMP- 


AXIS- 

BATHY- 


BURNS- 
CALMA  485- 


CDC- 

DATE- 


DAWHAT- 

DKLOCATE- 

DKOPEN- 

DKREAD- 


DKWRITE- 


GCDIST- 


Punction  used  in  MINCON,  MERFIX  and 
RHUMB  to  calculate  meridional  parts 
for  the  latitude  component. 

Calcomp  plotter  subroutine  to 
automatically  scale  and  draw  axes. 

Subroutine  v;hich  determines  which 
quadrant  the  rhumb  ?.ine  will  pass 
through,  extracts  the  gridded  data 
and  calculates  the  profile  for  each 
MSQLOC  area . 

See  SYNC0N2R 

(1)  A  large  bed,  graphical  analog 
digitizer  manufactured  by  the  CALMA 
Corporation. 

(2)  A  processor  program  for  (1)  that 
initially  scales  the  synthetic  track 
from  charts. 

Control  Data  Corporation 

COMPASS  off-line  subroutine  which 
automatically  calculates  an  integer 
day,  month,  year  from  the  computer's 
interval  clock. 

See  SYNCHEX 

Subroutine  which  positions  read/write 
head  at  specified  relative  address. 

Subroutine  which  opens  disk  file. 

Subroutine  which  reads  blocks  of  data 
from  the  disk  file  in  groups  of  32 
words  or  larger. 

Subroutine  which  writes  blocks  of 
data  on  to  disk  file  in  groups  of  32 
words  or  larger. 

Subroutine  used  by  SEAARCH  for  direct 
solution  of  the  great  circle. 
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Precedias  page  blank 


GCPATH- 


Subroutine  used  by  SEAARCH  for  indirect 
solution  of  the  great  circle. 


GEODIST- 

GEOPATIi- 

GRIDBLK- 

LOOKUP- 

LINE- 

MERFIX- 

MINCON- 

MSQFQ- 

MSQLOC- 

PLOT- 

PLOTS- 

PUNOUT- 


Subroutine  for  the  direct  solution  of 
the  geodesic. 

Subroutine  for  the  indirect  solution 
of  the  geodesic. 

Subroutine  which  calls  in  the  gridded 
data  from  the  random  access  storage 
device  for  BATHY, 

Subroutine  which  "looks  up"  or  extracts 
the  relative  address,  block  size  and 
the  column  and  row  information  for  each 
MSQLOC  area  from  the  random  access 
storage  device  previous  to  passing  this 
information  to  BATHY. 

Calcomp  plotter  subroutine  to  automati¬ 
cally  draw  a  line  as  a  function  of 
X  and  y. 

Subroutine  which  calculates  the  rhumb 
line  distance  and  sets  up  a  scaling 
factor  for  nautical  miles  along  a 
profile. 

Subroutine  used  to  calculate  the  start 
point  for  a  profile  within  a  MSQLOC 
area. 

Subroutine  used  to  calculate  in  part 
the  MSQLOC  area  numbers  for  points  on 
the  profile  path. 

Marsden  Square  Locator  Number  (Marsden 
square  system  is  a  numbered,  10  degree 
rectangular  grid  of  the  world  which  is 
subdivided  further  into  5  and  1  degree 
squares) . 

Calcomp  plotter  subroutine  which  moves 
pen  in  x  and  y  direction. 

Calcomp  plotter  subroutine  which 
initiates  plotter  action. 

Subroutine  which  places  each  MSQLOC 
area  profile  on  magnetic  tape. 
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RHUMB- 


SEAARCH- 

SPLico::^- 

SPLINE- 

SPLINT- 

SYMBOL- 

SYNBAPS- 

SYNBAPSl- 

SYNBLOCK" 

SYNCARD" 

SYMCHEX- 

SYNCOK2R- 

SYNGRID- 

SYNPLOT- 


Subroutine  using  AMP  to  compute  the 
rhumb  line  (approximation  of  a  chord 
of  a  great  circle  on  a  Mercator  pro¬ 
jection)  bearing  through  an  MSQLOC 
area. 

Subroutine  used  to  generate  a  great- 
circle  path. 

Subroutine  used  by  SPLINE  for  cubic 
spline  calculations. 

Subroutine  for  the  cubic  spline 
aJgorithro. 

See  SYNGRID 

Calcomp  plotter  subroutin.o  which  plots 
alphanumeric  characters  and  symbols. 

Synthetic  Bathymetric  Profiling  System. 

Accessing  program  which  produces  a 
depth  range  profile  on  i.iagnetic  tape 
for  each  MSQLOC  area . 

Program  v/Viich  loads  griclded  bathymetric 
data  into  the  random  access  storage 
device. 

Program  which  checks  longitude  of  data 
points  and  depth  values. 

Program  which  track  plots  data  points 
on  a  Mercator  projection  at  the  scale 
of  the  source  manuscript. 

Program  which  plots  contours  of  the 
gridded  data  on  a  Mercator  projection 
at  the  scale  of  the  source  manuscript. 

Program  v;hich  transforms  synthetic 
track  line  data  into  gridded  bathy¬ 
metric  date!  at  sev'en  points  per  card. 
This  is  the  primary  structuring 
program - 

Accessing  program  V7hich  links  together 
the  profiles  on  magnetic  tape  produced 
by  SYNBAPSi  for  each  M.SQLOC  area  to  plot 
a  great  circle  profile.  This  program  is 
usually  run  linked  to  SYT^BAPSl. 


SYNSTAT- 


SYNTABLE- 


SYNTRACK- 


TIMELEFT- 


UTILITY- 


Status  program  which  queries  random 
access  storage  device  for  listing  of 
file  key,  relative  address,  block 
size,  number  of  rows  and  columns  and 
date  that  data  were  added  to  storage 
and/or  actual  gridded  data. 

Traffic  director  program  which  supplies 
relative  address,  block  size  and  file 
key  to  SYNBLOCK  for  the  accurate  place¬ 
ment  of  blocks  of  gridded  bathymetric 
data  on  the  random  access  storage  device. 

Program  which  outputs  header,  track, 
data  and  blank  cards  and  conducts 
error  checks.  Input  is  a  scaled  data 
tape  from  the  CALMA  485-  processor 
program. 

COMPASS  on-line  function  which  extends 
time  mark  from  computer's  interval 
clock. 

Systems  program  which  loads  gridded 
bathymetric  data  cards  on  mangetic 
tape. 
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APPENDIX  A 


Preparation  of  Charts  for  Digitization 


The  5-degree  square  unit,  around  which  the  data  base  is 
created,  has  been  explained  in  the  "Outline  of  the  System"  and 
in  figure  3,  4,  5,  and  6.  Paper  copies  of  the  contour  charts, 
which  are  on  a  Mercator  projection,  are  used  to  prepare  the 
basic  manuscripts  for  digitizing.  Sufficient  overlap  around 
each  5-degree  square  is  required  to  provide  5  minutes  on  all 
sides  for  the  MSQLOC  area  and  an  additional  5  minutes  on'  all 
sides  for  interpolation  of  the  track  input  data  (fig.  6) ;  The 
manuscript  size  is  then  at  Teast  320  minutes  by  320  minutes 
regardless  of  the  chart  scale.  Ideally,  the  nanuscript  slioulc' 
consist  of  one  easy-to-handle  document.  However,  because  chart 
formats  vary,  this  is  not  alv/ays  oossible.  A  case  in  point  is 
the  addition  of  large  scale  survey  of  a  nev.’ly  discovered  seamount, 
to  a  regional  chart. 

One  method  of  handling  this  is  to  digitize  the  two  charts 
separately,  then,  substitute  the  synthetic  tracks  from  the  new 
seamount  chart  for  those  in  the  corresponding  section  of  the 
older  regional  chart.  A  second  method  is  to  prepare  a  contour 
selection  overlay  for  the  seamount  chart,  photographically  reduce 
it  to  the  scale  of  the  regional  chart,  make  a  print  at  that  scale, 
attach  the  print  to  the  regional  chart  and  match  the  contours. 

This  method  also  can  be  used  with  transparent  media. 

The  smallest  cell  selected  for  SYNBAPS  is  a  5-minute  (meridional 
part)  square  with  a  depth  value  at  the  four  corner  intersections. 

The  synthetic  tracks  of  input  depth  points  are  usually  taken  at  a 
5-minute  spacing  on  a  Mercator  projection.  In  high  frequency  data 
areas,  additional  tracks  of  data  at  1-,  2-,  3-,  or  4-minute  spacing 
can  be  input  so  as  to  improve  the  four  cell  depth  values.  However, 
there  is  a  limit  to  how  much  improvement  can  be  made  without 
losing  some  of  the  high-frequency  detail.  One  improvement  would 
use  a  smaller  cell  size,  but  this  makes  random- access  storage  device 
data  storage  requirements  very  large.  Thus,  small  features  that 
fall  within  a  5-degree  cell  can  be  lost  to  the  data  base, 
especially  if  they  are  not  picked  up  at  the  input  or  structuring 
phase. 

It  is  necessary  to  interpolate  the  beginning  and  end  points 
for  each  track  in  the  overlap  areas.  This  is  not  a  requirement 
for  short  tracks  within  the  body  of  the  MSQLOC  area.  These 
points  may  be  visually  interpolated  by  the  analyst  or  by  an 
experienced  digitizer  operator.  This  interpolation  need  only 
be  to  the  nearest  20  fathoms  or  about  one- tenth  the  contour 
interval. 


A-1 


The  output  from  the  syNC0N2R  program  is  a  contour  plot  of 
the  MSQLOC  area.  Although  this  output  is  not  a  primary  product 
of  the  system,  it  is  used  for  checking  and  may  be  a  useful 
byproduct  as  rough  automated  contours.  Because  of  the  S-minute 
cell  size  and  the  nature  of  the  interpolation  scheme,  large 
flat  areas  tend  to  break  up  on  the  contour  plot.  This  break  up 
of  contours  is  not  an  error  in  the  data  and  does  not  affect 
the  profile  generation.  To  improve  the  contour  output  aesthe¬ 
tically,  the  interpolation  can  be  improved  by  adding  contours 
in  key  locations.  In  areas  of  rough  topography  this  improvement 
will  not  be  necessary.  The  first  example,  around  seamounts  or 
a  seamount  group,  is  shown  in  figure  A-1.  Usually  the  added 
contour  is  placed  outside  the  base  contour  to  cutoff  or  terminate 
the  interpolation  adjoining  a  flc.t  area  or  to  define  the  sea¬ 
mount  base.  The  second  example  is  for  domes,  rises,  ridges  or 
tablemounts  (fig.  A-2) .  Here  the  added  contours  are  on  the  top 
of  the  structure  in  order  to  cutoff  or  terminate  the  interpolation 
on  their  flat  or  gently  rounded  summits.  The  third  example  is 
for  noses  or  spurs  (fig.  A-3) .  Although  this  feature  is  similar 
to  those  in  figures  a-1  and  A-2,  short  disconnected  contours  may 
be  needed  if  the  spur  slopes  are  gentle.  In  all  these  examples, 
the  track  direction  was  assumed  to  be  left  to  right. 

The  boundary  condition  is  a  special  case  of  endpoint  inter- 
pdlation.  Whenever  an  island  or  continent  is  encountered,  the 
zero  contour  or  sea  level  is  handled  as  shown  in  figure  A-4. 

On  the  SYNC0N2R  program  the  zero-contour  level  should  never  be 
plotted,  but  the  1-fathom  or  1-meter  contour  should  be  inter¬ 
polated  to  show  the  coast  line.  In  profiling,  the  punched  card 
depth  values  after  the  first  zero  usually  are  discarded  and  the 
profile  terminated  at  that  range. 
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FIGURE  A-2.  ADDED  CONTOURS  ON  DOMES  OR  RISES 
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FIGURE  A-3.  ADDED  CONTOURS  AROUND  A  SPUR 
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FIGURE  A-4.  BOUNDARY  CONDITIONS  FOR  ZERO  CONTOUR 
LFv/EL 
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APPENDIX  B 
FORTRAN  Programs  for 
Structuring  SYNBAPS 


All  programs  and  subroutines  listed  in  this  appendix  are 
subject  to  change  without  notice.  Modifications  within  the 
programs  and  adoption  of  the  system  for  other  computers  will 
necessitate  major  changes.  The  author  should  be  contacted  for 
the  most  recent  versions  oc  these  program  .  . 
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APPENDIX  C 
FORTRAN  Programs  for 
Accessing  SYNBAPS 


All  programs  and  subroutines  listed  in  this  appendix  are 
subject  to  change  without  notice.  Modifications  within  the 
prcgic.ms  and  adoption  of  the  system  for  other  computers  V7ill 
necessitate  major  changes.  The  author  should  be  contacted  for 
the  most  recent  versions  of  these  programs. 
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MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST55I0.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 .  In  accordance  with  reference  (a),  a  declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A;  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 
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